Statistics of resonance width shifts as a signature of eigenfunction nonorthogonality 
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We consider an open (scattering) quantum system under the action of a perturbation of its closed counterpart. 
It is demonstrated that the resulting shift of resonance widths is a sensitive indicator of the nonorthogonality of 
resonance wavef unctions, being zero only if those were orthogonal. Focusing further on chaotic systems, we 
employ random matrix theory to introduce a new type of parametric statistics in open systems, and derive the 
distribution of the resonance width shifts in the regime of weak coupling to the continuum. 
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The classical question of how energy levels of a quantum 
system get shifted under the action of a perturbation kept at- 
tracting renewed attention during the last two decades, mostly 
due to the established universality of such a parametric mo- 
tion for systems with chaotic dynamics or intrinsic disorder 
iHHI]- In particular, the distributions and correlation func- 
tions ofparametric derivatives of energy levels ("level veloc- 
ities") tl|-r3[] and their second derivatives ("level curvatures") 
J3l were found explicitly using the methods of random matrix 
theory (RMT) fl, and also verified, e.g., in microwave bil- 
liard experiments Jg]. The other reason for such an interest is 
the recent development of the fidelity concept as the measure 
of the susceptibility of internal dynamics to perturbations 0]. 

Experimentally, the energy levels are mostly accessible by 
means of a scattering setup [8J. From such a viewpoint, para- 
metric dependencies of scattering characteristics, like phase 
shifts and time delays Jit], conductances iflf^l and S matrix 
elements 111 111 were under intensive study. As to the parental 
discrete energy levels, they are converted into the resonances 
with finite lifetimes, since the original closed system becomes 
open (unstable). Such resonances manifest themselves in the 
energy-dependent S matrix as its poles in the complex energy 
plane, and can be analytically described as the complex eigen- 
values of an effective non-Hermitian Hamiltonian Il2l - ll4ll . 
Notably, the corresponding eigenfunctions are not orthogonal 
in the conventional sense but rather form a biorthogonal sys- 
tem. Their nonorthogonality plays an important role in many 
applications, e.g., describing interference in neutral kaon sys- 
tems |15p. influencing branching ratios of nuclear cross sec- 
tion s Ml 611 . and yielding excess noise in open laser resonators 
[17l[l8|]. It also features in decay laws of quantum chaotic 
systems lfl9ll and in dissipative quantum chaotic maps lEoll . 

In such a context the question of parametric motion of reso- 
nances and associated resonance states in open systems arises 
very naturally, but to the best of our knowledge has never 
been properly addressed. Our goal here is to begin filling in 
that gap by considering universal statistics of the shifts in the 
resonance widths under a generic perturbation in chaotic sys- 
tems. In particular, we will demonstrate that such shifts are a 
clear manifestation of eigenstate nonorthogonality, thus pro- 
viding a promising way to probe this spatial characteristics by 
purely spectroscopic tools. To this end let us stress that statis- 



tics of complex poles and lifetimes in chaotic systems can be 
verified via accurate scattering experiments in microwave bil- 
liards 12111 or photonic crystals |22(]. It can be also extracted 
from realistic computer simulations of quantum graphs IT23ll . 
semiconductor superlattices 12411 . dielectric microresonators 
[25] or system of randomly interacting fermions IT26Tl . As to 
the theoretical framework, it mainly relies on studying the rel- 
evant non-Hermitian RMT IU3[ Il4ll . understanding of which 
has substantially improved over the last two decades; see, e.g., 
Ref. 112711 and references therein. Note that the spatial proper- 
ties related to the associated bi-orthogonal eigenvectors are 
known to a much lesser extent (l8 , 28 - 30l] ■ 



As is well known H 121 - 41411 (see also 113 111 for most recent 
reviews), resonance phenomena involving a group of N inter- 
fering resonances can be adequately described in terms of the 
effective non-Hermitian Hamiltonian 



H eS = H- i-fWW^ , 



(1) 



which governs the dynamics of the open system. Here, the 
Hermitian NxN matrix H corresponds to the Hamiltonian of 
the closed counterpart, whereas the entries W£ of the rectan- 
gular NxM matrix W are the decay amplitudes that describe 
coupling of N discrete energy levels, labeled by n, to M de- 
cay channels, labeled by c. The coupling strength to the con- 
tinuum is controlled by the dimensionless positive constant 7, 
with 7 <C 1 (7 = 1) being the particular case of weak (perfect) 
coupling. The eigenvalue problem for the full non-Hermitian 
matrix H c h reads as follows 



Mcs\Rn) = £n\Rn) and (L n \H e {£ = £ n {L r 



(2) 



and determines its complex spectrum £ n = E n — §r„. Here 
E n stands for the resonance positions (energies) and r„ > 
denotes the corresponding widths |32ll . The set of the asso- 
ciated left and right eigenvectors (resonance wavefunctions) 
satisfies the conditions of biorthogonality, (L n \R m ) = 8 nm , 
and completeness, XL =1 \R n ){L n \ = 1. 

In such a framework a perturbation of the closed counter- 
part can be modeled by the term aV, where V is a Hermitian 
NxN matrix and the real constant a is to control the pertur- 
bation strength. The resonance parameters for the perturbed 
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open system are then to be determined from solving the fol- 
lowing spectral problem for the right eigenstates, 



{% e n + aV)\B! n )=E' n \B! n 



(3) 



and a similar problem for the left eigenstates (L' n \, In the 
case of the weak perturbation, a -C 1, one can follow the 
standard perturbation theory routine and seek each of the two 
eigenvectors corresponding to the new eigenvalue £' n as an 
expansion over non-perturbed basis of H^ff, with necessary 
modifications induced by bi-orthogonality 113 3Q . To the first 
order in a this readily yields the expression for the shift of the 
?ith resonance in the form 



S£ n =£' n — £ n = a{Ln\V\Rn 



(4) 



generalizing the standard result to the non-Hermitian case. 

The resonance shift contains both real and imaginary 
parts, since (L n \ ^ (|i?„))^ in general 13411 . At this point we 
stress that a nonzero value of the imaginary part of S£ n is in- 
duced solely due to nonorthogonality of the resonance states. 
This fact becomes apparent in the following equivalent repre- 
sentation for the resonance width shift ST n = —2lm(S£ n ): 



ST n = ia((L n \V\R 



(Rn\V\L n )) 



(5) 



where G n stands for the following Hermitian operator: 

G n = WW f \n){n\V + V\ri)(n\WW ] . (7) 

Aiming to describe the universal statistics of the resonance 
shifts for generic chaotic systems (e.g., open billiards, quan- 
tum dots, quantum graphs), we follow the standard paradigm 
Jil Ht] and model H by a random N x N matrix drawn from 
the Gaussian orthogonal (GOE) or unitary (GUE) ensemble, 
depending on the presence or absence of time-reversal sym- 
metry, respectively. Universal fluctuations are then expected 
to occur in the limit N 3> 1 at the local scale of the order of 
the mean level spacing A ~ 1/N. Without loss of general- 
ity, one can restrict the consideration to the spectrum center, 
where A = Xtt/N and 2A is the semicircle radius (A needs 
to be rescaled if E 7^ 0). As to the coupling amplitudes, they 
can be taken 0] as real (GOE, /3=1) or complex (GUE, 0=2) 
Gaussian random variables with zero mean and the variance 
( w n w m) = (A/7r)<5 m „<5 Q \ the final results being model- 
independent provided the number of channels M -C N. This 
readily yields the well-known result for the resonance width 
distribution in terms of \ 2 V distribution with v = Mf3 degrees 
of freedom (Porter- Thomas expression at M = 1 and (3 = 1), 



T(M0/2) 



M0/2-1 -0k/2 



(8) 



where we have made use of the completeness condition to ex- 
pand \L n ) = Y^rn \ R m) (L m \L n ) and also introduced V nm = 
(R n \V\R m ) and U nm = (L n \L m ). Since by construction 
V nm = y* nn and U nm = U* m , only the terms with to ^ n 
contribute to the sum (0. The matrix U is just the Bell- 
Steinberger nonorthogonality matrix IU5I1 (see also a compact 
description in IU3I0 . and U nm ^ S nm in general. Thus, the 
resonance width shift ([5]l would generically vanish only if the 
resonance states were orthogonal ll35ll . being thus a sensitive 
indicator of their nonorthogonality. 

In the rest of the Letter, we concentrate on the regime of 
weak coupling to the continuum, 7 <C 1, which permits 
complete analytical investigation, and is the one that is the 
most easily realized experimentally. Under such an assump- 
tion the non-Hermitian part i^WW^ of H c s can be treated 
as the perturbation of the Hermitian part H . To the leading 
order in 7 the resonance positions coincide with the energy 
levels of the closed system, H\n) = E n \n), whereas the reso- 
nance widths are given by expression T n = 2^ (nlWW^ \n) = 
27 J2c=i l^nl 2 - Similarly, the right eigenvectors of H e e are 
readily found to be \R n ) = \n) - «7E m /„ ^T-eZ" 1"^' 
while the left ones read (L n 
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Substituting such expressions into Eq. ©, one finds that to the 
leading order in both a and 7 the energy shift is given by the 
standard expression for the closed system SE n = a(n\V\n), 
whereas the shift in the resonance width is determined by 



ST n = 2« 7 E HG "' m) 



E n — E TI 



(6) 



where n n = |^ stands for the width measured in units of 
the mean partial width (per channel). Distribution $Q has the 



mean value 



M and variance var(«) 
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In the limit N 3> 1, the rescaled matrix elements i> 
N(m\V \n) I y/Tv(y^) of the perturbation in the eigenbasis of 
H become normally distributed random variables (13611 . This 
results in the Gaussian distribution of the energy shifts 5E n 
(i.e. "level velocities") with var(£E„) = j^Tr(V 2 ). At the 
same time the width shifts 8T n must clearly have much less 
obvious distribution due to the nontrivial structure of Eq. ([6j. 
To find the latter distribution explicitly, we first scale the width 
shifts in the natural units to get rid of the model-dependent 
features, and introduce 



5T r , 



tin 



2iy/2Pvar(6E n ) 



(9) 



In close analogy with the case of the closed systems such a 
rescaling is expected to capture universal (local) fluctuations 
related to the parametric motion in generic open systems. At 
the next step it is instructive to follow Ref. 13011 and treat the 



scalar product A J K (WW^) mn — Zm as a projection of the 

^/-dimensional vector of the decay amplitudes ^p£W m onto 
the direction determined by the vector W n at given n 7^ to. 
As a result, the expression for the rescaled width shifts takes 
the following convenient form: 
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The advantage of such a parametrization is that the corre- 
sponding modula and angles turn out to be statistically in- 



(n) 

dependent [13]. The projections z m can then be shown to 
be normally distributed random variables (real at /3 = 1 and 
complex at j3 = 2) at any M > 1 l30ll. 

We now compute the probability distribution Vm (y) of the 
rescaled width shifts (at the spectrum centre) defined as 



JY 



V M (y) = A{J2S(E n )6(y 



(11) 



where the angular brackets denote the spectral average over 
both energies E n and widths K n , whereas the bar stands for 
the averaging over the quantities Zm and utl all of them 
being statistically independent. The latter task can be easily 
performed by considering the Fourier transform, yielding 



n 



\E n — En 



[(£„ - £„J 2 + ^(wA/ttV?) 2 ]* 3 / 2 



(12) 



Substituting this expression back to Eq. (fTTT i. one can then 
make use of the known explicit form of the joint probability 
function of all eigenvalues to integrate out the nth eigenvalue. 
The distribution can be finally represented as follows 



V M {y) = 



(13) 



where the width distribution P K; ' (k) is given by Eq. ([8]) and 
the function ^(y) = ^e^C^l^uj) is the Fourier- 
transform of the following ratio of spectral determinants for 
the (GOE or GUE) matrix Hi of the reduced size N — 1: 



N-l 



det(#i) 2 ' 3 



det[H 2 + (uA/^)f /; 



(14) 



The constant c^'_ x = (det \Hi\~P) ensures Cj^ii(O) = 1 
that automatically guarantees the normalization of distribution 
(fT~3T > to unity. Objects similar to Eq. (fT4l naturally arise in the 
analysis of weakly open chaotic systems as a result of separat- 
ing independent fluctuations in spectra and in wavefunctions. 
Following the methods developed in IU8I1 for /3 = 1 and in 
1 3711 for P = 2, we have been able to calculate the limiting 
expressions at N S> 1, with the explicit result being 



c^M = i [MA^(M) + w 2 a- 2 (M 



c£\oj) = e~ 2 M ( 1 + 2M + ±f + 



(15) 



where K v [x) stands for the modified Bessel (Macdonald) 
function. Taking the Fourier transform, we finally arrive at 



y 



(y) = 



6(1 

35 H 



hy 2)5/2 ' 

Uy 2 + 3y 4 



127r(l + y 2 ) 4 



0=1 

P=2 



(16) 
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FIG. 1 . Distributions of the resonance width shifts for weakly open 
chaotic systems with preserved (/3=1, top) or broken (/3=2, bottom) 
time-reversal symmetry at M — 1 (•), 2 (o), 5 (★) and 10 (□) open 
channels. The solid lines show our analytical prediction, Eqs. \13\ 
and J16K The symbols stand for numerics with 2000 realizations 
of 250x250 GOE (,3=1) or GUE (fi=2) random matrices (only 25 
levels around the spectrum centre were taken for each realization). 



Combination of Eq. dT~3T > with Eqs. (O and ( [ToT i completely 
solves the problem of universal statistics of the resonance 
width shifts in a chaotic system weakly coupled to the con- 
tinuum via M equivalent channels. We see that far tails of the 
distribution decay as Vm(v) oc y _ (' 3 + 2 ) due to the Wigner- 
Dyson level repulsion at small energy level separations, the 
feature which such a distribution shares with that for level 
curvatures J3]. In contrast to the level curvature distribution, 
the broadness of the width shift distribution dTJb can be ad- 
ditionally controlled and is proportional to \fM ~ yjvai{n). 
Physically, this gives the variance of widths the role of a uni- 
versal parameter that controls the degree of nonorthogonality 
in weakly open chaotic systems |2^, 3^1 . In the limit of many 
weakly open channels M ^> 1 (but still M <C N) the widths 
cease to fluctuate, so distribution ^ becomes very narrow and 
peaked around its mean value (k) = M. In such a limit the 
width shifts are still widely distributed, with the probability 
density for the scaled variable?; = y/y/M being givenjustby 
the function <f>^> (y), Eq. (TToT l. These findings are illustrated 
on Fig. 1 which shows the distribution Vmiy) at several val- 
ues of M. Also shown are the results of straightforward nu- 
merical simulations of the width shifts (O with GOE/GUE 
random matrices, the agreement being flawless. 
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In summary, we have shown that the change in the reso- 
nance widths due to external perturbation is a sensitive in- 
dicator of eigenfunction nonorthogonality in open quantum 
or wave systems, and have computed its distribution analyti- 
cally for weakly open chaotic systems with preserved or bro- 
ken time-reversal symmetry. Our predictions can be veri- 
fied, e.g., in experiments with reverberant dissipative bodies 
Jilt], microwave cavities ll39ll . and elastic plates JJjJ, where 
other aspects of nonorthogonality were investigated. Poten- 
tially, they may also have implications for the other type of 
non-Hermitian systems, those with PT-symmetry, which is a 
rapidly developing field |41fl . We also note a link to a more 
formal concept of pseudospectra of non-selfadjoint operators 
considered mostly in mathematical literature [1421] . The gener- 
alization of our results to the case of arbitrary modal overlap 
is another challenging problem to consider in future studies. 
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